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lateral cyclic pitch, deg 

longitudinal cyclic pitch, deg 
blade section drag coefflcient 
blade section lift coefficient 
blade section moment coefficient 
QIpiORyitR^ 

TlpiWi)^irR^ 

PFIp iSlRYnR^ 

rotor equivalent parasite drag area, m’ 
lift force 

moment; Mach number 

moment about blade feathering axis, N-m 

rotor roll moment, N-m 

rotor pitch moment, N-m 

rotor power, hp 

rotor propulsive force, N 

rotor torque, N-m 

rotor blade radial station (r = 0 at center of rotation, r = R at tip) 
rotor radius 

rotor thrust, N pn rv> , 
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rotor side force, N 

blade section angle of attack, deg 

rotor coning angle, deg 


NOT 
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cyclic flap, pitch (forward) of tip path plane, deg 
Pis cyclic flap, roll (to left) of tip path plane, deg 

mean blade lag angle, deg 
elastic torsion 

fl .7 5 rotor collective pitch angle (at r/Jt = 0.75), deg 

Xypp rotor inflow velocity (divided by tip sp'^ed), in the tip path plane 

fi rotor advance ratio: helicopter forward speed divided by rotor tip speed 

p air density 

a rotor solidity: ratio of total blade area to rotor disk area 

^ rotor blade azimuth angle, measured from downstream direction 

rotor rotational speed 
( ) time derivative 




COMPARISON OF THREE METHODS FOR CALCULATION OF HELICOPTER 


ROTOR BLADE LOADfWG AND STRESSES DUE TO STALL 


Wayne Johnson 

Ames Research Center 
and 
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SUMMARY 


A comparison is made cf the results of three methods for calculating the effects of dynamic 
stall on the performance, airloads, and blade stresses of a helicopter rotor at high loading. The three 
dynamic-stall methods considered predict essentially the same performance and trim for the rotor. 
TTiey give roughly the same mean bending moments, but the peak-to-peak torsion and bending 
moments differ by 25 to 40 percent, and there are differences in the details of the predicted blade 
motion and stresses. The latter are due to significant differences in the dynamic stall aerodynamic 
loads, particularly the aerodynamic pitch moment, predicted by the three methods. 


INTRODUCTION 


Rotor blade stall usually is the limiting factor in helicopter maneuver and high-speed capabil- 
ity. Particularly important are the vibration and control system loads due to the high transient lift 
force and pitching moment of the periodically stalling blade. An adequate theoretical model of the 
stall process is therefore required for accurate design of the helicopter and prediction of its perfor- 
mance and limitations. Prediction of the aerodynamic loads of the stalled rotor blade is com- 
plicated, however, by the complex aerodynamic environment in which the blade operates. Rotor 
blade stall is always an unsteady and three-dimensional phenomenon, so it is not sufficient to use 
static, two-dimensional airfoil section data. 

Considerable theoretical and experimental research has been coi.ducted into the nature of 
unsteady, or dynamic, stall. Based on this research, several semi-empirical methods have been 
developed to incorporate the effects of dynamic stall into the calculation of the performance, 
airloads, and blade stresses of the helicopter rotor at high loading. 

This report compares the results of three such methods, which are described briefly below. 
Additional details are provided in appendix A and the references cited. 

1. UARL Method (refs. 1 ,2). Tabular data for cg and as d function of a, a, and a are used 
to calculate the stall loads. These data were synthesized from oscillating airfoil test results. Static 
airfoil data are used below the stall angle of attack. 



2. MIT Method (refs. 3-S). High transient cg and nosedown are used at the occurrence of 
dynamic stall. The dynamic stall loading is basically impulsive, and then static stall values are used. 
The peak cg and c^, as a function of a, are based on the experimental results of reference 3. Static 
airfoil data are used below the stall angle of attack. 

3. Boeing Method (refs. 6,7). The actual angle of attack is corrected as a function of d to 

obtain an effective angle ofjyn- Then cg, f j, and are obtained from static airfoil data using a^yn. 

The angle-of-attack correction is based on oscillating airfoil test results. 


CALCULATION CASES 


The stall methods were compared by calculating the rotor performance and blade loads at a 
high speed, high thrust (jx = 0.333, Cjla = 0.09) operating condition. Appendix B gives the details 
of the rotor operating state; the geometrical, inertial, stnictural, and aerodynamic characteristics of 
the rotor blade; and the parameters involved in the calculation process. These rotor and operating 
condition parameters were developed by Ormiston (ref. 8), for a general comparison of several 
methods of predicting the loads on a hypothetical rotor. The rotor performance and loads were 
calculated using a computer program based on reference 9, and incorporating the three stall 
methods considered here. Hence, the only difference in the methods is in the calculation of cg, 
and for the blade sections, and this difference is apparent mainly at high angle of attack since all 
three methods use the same static airfoil data for low angle of attack. The following cases were 
considered: 

CASE 1: Linear, incompressible aerodynamics (cg = 5.7«, c^ = 0.0086, = 0); in other 

words, no stall. This base case is also used to examine the influence of nonuniform 
inflow (Case lA, uniform inflow) and the influence of elastic blade motion 
(Case IB, which uses only rigid flap and lag modes, as well as uniform inflow). 

CASE 2: Static stall; static airfoil data used for all a 

CASE 3: UARL stall method. 

CASE 4. MIT stall method. 

CASE 5: Boeing stall method. 

Static airfoil data for an NACA 0012 section, in tabular form, wereused in all three methods for the 
loads below stall. All the cases (except lA and IB) used nonuniform inflow, which was calculated 
by a separate program. Hence, the inflow distribution was prescribed in these calculations, and was 
the same for every case. 


RESULTS AND DISCUSSION 


Figure 1 shows the nonuniform inflow distribution used in the calculations. The role of the 
inflow distribution is indicated in Figures 2 and 3, which are polar plots of the angle-of-attack 
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distribul.v .1 over the rotor disk for cases lA and 1, respectively. Case I A uses uniform inHow and 
case I the nonuniform inflow; neither includes stall aerodynamics or the effect of the blade motion 
due to stall. 

Table I presents tlie performance and trim data calculated for the seven cases considered. 
Cases I, I A, and IB are with no stall; cases I A and IB examine the influence of uniform inflow and 
rigid blade motion, respectively. Case 2 uses static stall. Compared with the no-stall and the 
dynamic-stall cases, the static stall case shows very high values for power required, rotor flapping, 
and control angles (collective and cyclic), even though the thrust is still somewhat below the target 
value. These results imply that this operating condition (the combination of thrust, propulsive 
force, and forward speed) would be beyond the capability of this rotor if the blade stalled statically. 
The loads and the bending and torsion moments calculated for case 2 are unrealistic, therefore, and 
are not included in the other comparisons presented here. The three dynamic stall methods (cases 3, 
4, and 5) give essentially the same results for performance and trim. Comparing, for example, the 
power required for these three cases with that for the no-stall case, it is evident that the stall effects 
are quite important. However, the effects predicted using the dynamic-stall models are much less 
than those obtained using static-stall data (case 2). 

Table I also presents the results for the vertical and inplane root shear forces of the rotor. 
Only the harmonics transmitted to the helicopter body are given. There is a trend to increased 
magnitude of the higher harmonics due to stall. The variation of the root shear forces between the 
cases is an indication of the general sensitivity of the problem of calculating helicopter rotor 
vibration, as well as the difference between the stall models. 

Figure 4 compares the elastic torsion motion of the blade, for the three dynamic stall methods. 
Since only one mode is used to represent the torsion motion, the torsion moments are directly 
proportional to the torsion deflection. Figure 4 shows significant differences in the details of the 
torsion moment waveform, particularly on the retreating side where dynamic stall is involved. The 
peak-to-peak torsion moments given by the three methods differ by about 40 percent. Figures 5 
through 8 show the peak-to-peak and mean flatwise and edgewise blade bending moments for the 
three dynamic stall cases and the no-stall results (case 1). The mean bending moments predicted by 
the three stall cases are roughly the same. The higher mean edgewise moment compared wi^h the 
no-stall case is due to the increased blade section drag associated with stall. Consequently the 
increase in the edgewise bending moment is about the same as the increase in the power required. 
The calculated peak-to-peak bending moments show roughly the same general characteristics, but 
the magnitude predicted by the three stall methods differ by about 25 to 40 percent. 

The differences in the calculated blade torsion and bending moments may be traced to the 
blade aerodynamic loading. Figures 9 to 12 present the section lift force, and figures 13 to 16 the 
section pitching moment for the no-stall case (case 1) and the three dynamic-stall methods (cases 3, 
4, and S). There are significant differences in the aerodynamic loading, particularly in the pitch 
moment, due to the fundamental differences in the calculation of the dynamic stall loads in the 
three methods. 



CONCLUDING REMARKS 


The three dynamic-stall methods considered predict essentially the same performance and trim 
for the rotor. They give roughly the same mean bending moments, but the peak-to-peak torsion and 
bending moments differ by 2S to 40 percent, and there are differences in the details of the pre- 
dicted blade motion and stresses. The latter are due to significant differences in the dynamic-stall 
aerodynamic loads, particularly the aerodynamic pitch moment, predicted by the three methods. 
The semi-empirical methods that have been developed to incorporate dynamic-stall effects into 
helicopter airloads calculations are useful and important advances. However, the fundamental differ- 
ences in the aerodynamic loads predicted by the various methods suggest a need for further study 
into the fundamentals of the aerodynai.iic problem before a fully confident prediction of all 
helicopter rotor stall effects is possible. 

Ames Research Center, NASA 
and 

U.S. Army Air Mobility R&D Laboratory 
Moffett Field, Calif. 94035, July 10, 1974 
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APPENDIX A 


THREE METHODS OF CALCULATING DYNAMIC STALL EFFECTS 


UARL Method 

The UARL method calculating dynamic stall airloads is based on data from tests of a NACA 
0012 airfoil oscillating in pitch. The data for lift and moment coefficients cg and are correlated 
as a function of a, A =d(c/2F,and 5 = a(c/2K)’, by cross-t>lotting the oscillating airfoil data. The 
stall method uses these data for cg and as a function of a, A , and B. The data are given in tabular 
form in reference 2. Figure 1 7 presents typical results for cg and as a function of a, for several 
values of a and with of = 0. Compressibility effects are accounted for by scaling the incompressible 
data for cg and c^. 


MIT Stall Method 


The MIT stall method is a theoretical model of the dynamic-stall loading of an airfoil for the 
calculation of helicopter airloads and blade motion. It is based on the experimental data of ref- 
erence 3 for the peak transient airfoil loads during dynamic stall. When the angle of attack reaches 
the dynamic-stall angle (which is above the static-stall angle dynamic stall of the section 
occurs, resulting in a sharp increase in the lift and nosedown moment. The peak cg and are 
functions of ac/V at the instant of stall. The rise time of the loads to the peak values, and the fall 
time to static stall is small. Hence, the model involves essentially impulsive loading of the blade 

occurring at dynamic stall, with maximum loads eg and , after which the blade has 

msx *rndx 

static-stall loads until the flow reattaches (reattachment is assumed to occur at the static stall angle 
%)• 

When the blade section angle of attack exceeds dynamic stall occurs at that section, 
producing high transient loads. From reference 3, th? cg and at stall are: 


*max 


1 -I- 40l(ic/F| 
3 


6tclV< 0.05 
aclV> 0.05 


'■m 


max 


-0.15 

0.283- 21.7 lac/Kl 
,-0.80 


btclV< 0.02 
0.02 <ck/K<0.05 
6tclV> 0.05 


A linear rise of the lift and moment to the peak values, in time is used. Then i linear decay 

to static stall Cg and occurs. When the transient dynamic-stall loads have decayed, cg and 
are assumed to take the static-stall values. For the NACA 0012 airfoil, best correlation is obtained 
with«£)j= 15°, 10® anda^^ = 12® (for reattachment). 


Boeing Stall Method 


The Boeing stall method is a theoretical model for the unsteady aerodynamics of a blade 
section, including dynamic-stall and radial flow effects. It is based on static airfoil data (cg, and 
as function of a. and with the angle of attack corrected for unsteady aerodynamics: 


Cj = min 


^g^“dyn,L^ 



*equ 


~ '2 p 


whem A is the sweep angle, and is the skin friction drag coefficient. The correction for the 

effect of sweep on the lift is just (ignoring the angle-of-attack corrections for the moment) 
= cj^jfj/cos A. The lift coefficient is thus increased by the factor (cos AT‘ , with the restric- 
tion that the lift curve slope cannot exceed the unswept airfoil value The equivalent angle of 
attack otgq^ is the actual angle of attack corrected for unsteady potent!^ flow effects (see refs. 6 
and 7). The dynamic angle of attack is a dynamic-stall correction of a, due to the pitch rate 6c. 


“dyn = “ 


* yip s8n(a) 


The additional subscript I or M on a^y^ indicates whether the correction is for the lift or for the 

n.oment. The difference is in the constant A:, which for the NAC/V0012 airfoil section is {k in 
degrees): 


lift stall k = max (0, min (98, 1 24 - 1 6 IM)) 

moment stall k ^ max (0, min (64, 98 - 17 lAf )) 
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APPENDIX 0 


ROTOR CHARACTERISTICS. OPERATING CONDITIONS. 
AND CALCULATION PARAMETERS 


Rotor Physical Description 

Articulated blade with flap, lag, and feather hinges at hinge offset e; no ptecone, prelag, 
pitch/flap coupling, pitch/lag coupling, flap hinge damp^'r, or structural damping; spar centroid, 
aerodyanmic center, center-of-gravity, and elastic axis at quarter chord. 

Geometry 


radius, R 

7.62 m (25 ft) 

\ number of blades, N 

3 

i hinge offset, e/R 

0.04 

linear twist, from center of rotation to tip 

-10" 

chord, c 

0.558 m 

solidity, o 

0.070 

Inertial and structural parameters 


nmning mass, m 

11.97 kg/m 

blade weight, 

859 N 

first mass moment, 5^ 

320 k^m 

second mass moment, 

155*? kg-m* 

Lock number, 7 

8.43 

flatwise bending stiffness, 

8.61 XIO^ N*m* 

edgewise bending stiffness, 

1.148X10* N-m’ 

torsional moment of inertia, 

.004 kg-m* /m 

torsional strength, GJ 

8.61X10* N-m* 

Three flatwise modes (rigid flapping and two elastic bending modes), two edge 

and one elastic bending mode), and one toraon mode (elastic torsion, with a 

used in the calculation of the blade motion. 


Rotating natural frequencies of blade modes (per rev) 


Flatwiy* Edgewise 

Torsion 

1.030 0.24S 

5.836 

2.694 3.851 


; 5.055 


i 

1 Aerodynamic properties 


section airfoil 

NACA0012 

tip loss factor, B 

0.97 

root cutout 

0.15/; 


y 



Operating Conditions 


Rotor operating stages 
forward velocity, V 
tip speed, 

shaft angle-of-attack, 
advance ratio, /x 

advancing tip Mach number, 1 . 0,90 

Rotor force trim (target values) 
thrust, T 

propulsive force, FF 
side force, V 

helicopter parasite drag area, / 
rotor blade loading, Cj’/a 


148 knots 

229 m/sec (750 fps) 
0 

0.333 

0.895 


73,395 N (16,500 lb) 
8260 N (1857 lb) 

0 

2.32 m* (25 ft*) 
0.090 


Calculation Parameters 

Blade motion calculated by a harmonic analysis method, 

azimuth increment 5 ° 

number of harmonics 1 5 

tolerance on blade mode amplitude and 
velocity, in test for convergence 0.001 

number of blade segments (radial) 1 5 

Iterate on ^.7 s » 1 5 . and j to trim rotor thrust, propulsive force, and side fo 

tolerances; 

thrust 445 N 

propulsive force 44 N 

side force 44 N 
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TABLE I - PERFORMANCE AND TRIM DATA 
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TABLE 1.- PERFORMANCE AND TRIM DATA - Continued 
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TABL*^ I - PERFORMANCE AND TRIM DATA - Concluded 
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Figure 7.- One-half peak-to-peak edgewise blade bending 
moment. 



Figure 8.- Mean edgewise blade bending moment. 
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